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The details are presented of a new evolution algorithm for the characteristic initial-boundary value 
problem based upon an affine parameter rather than the areal radial coordinate used in the Bondi- 
Sachs formulation. The advantages over the Bondi-Sachs version are discussed, with particular 
emphasis on the application to the characteristic extraction of the gravitational waveform from 
C**) ' Cauchy simulations of general relativistic astrophysical systems. 



o 



X 



PACS numbers: 04.20.-q, 04.20.Cv, 04.20.Ex, 04.25.D- 



I. INTRODUCTION 



[ — I There has been important progress in computing accurate gravitational waveforms by means of Cauchy- 
characteristic extraction (CCE) [l|, whereby data from a Cauchy simulation provides the inner boundary data for a 
characteristic evolution extending to future null hffinty Z + , where the waveform is defined unambiguously. CCE has 
' ^j ' , become an important tool for gravitational wave data analysis Q. It has been applied to compute waveforms from 
simulations of binary black hole inspiral and mergers , from rotating stellar core collapse @, 0] , to explore the 
memory effect Q and to study the effect of spin on gravitational waves from precessing binary black holes Q . 
60- A CCE module [TO) has been prepared for public use as part of the Einstein toolkit lll|. The module is based 
upon the PITT null code [12, [lj|, which implements the worldtube-nullcone version [l4j of the Bondi-Sachs (lH 
|l6[ characteristic initial-boundary value problem. There are technical complications in applying the Bondi-Sachs 
formulation to CCE arising from the use of an areal radial coordinate to parameterize the outgoing null geodesies. 
^\ . This paper considers an alternative approach to the worldtube-nullcone problem which replaces the areal coordinate 
\& ' by an affine parameter. The details of an evolution algorithm for the affine system of Einstein equations arc presented. 
0^ , The comparative advantages with the Bondi-Sachs version for application to CCE are discussed. 

Recent success in simulating general relativistic astrophysical systems has been achieved by Cauchy codes, which 
(T) , evolve the spacetime metric inside an artificially constructed outer boundary. In doing so, it is common practice to 
compute the gravitational waveform from data on an extraction worldtube inside the outer boundary, using pertur- 
bative methods based upon introducing a Schwarzschild background in the exterior region. This has been carried 
out using the Regge- Wheeler- Zerilli [l?], [l8[ treatment of the perturbed metric, as reviewed in [19], and also by 
calculating the Newman- Penrose [13] Weyl curvature component ^4, as first done for the binary black hole problem 
in [2ll - |24j . In this approach, errors arise from the finite size of the extraction worldtube, from nonlinearities and from 
gauge ambiguities involved in the arbitrary introduction of a background metric. The gauge ambiguities might seem 
less severe in the case of "J 4 (vs metric) extraction, but there are still delicate problems associated with the choices 
of a preferred null tetrad and preferred worldlines along which to measure the waveform (see [25| for an analysis). 

In order to properly approximate the waveform at I + the extraction worldtube must be sufficiently large but at 
the same time causally and numerically isolated from errors propagating in from the outer boundary. Considerable 
improvement in the perturbative approach has resulted from techniques for dealing with large outer boundaries and 
extrapolating the extracted waveform to infinity. However, this is not an ideally efficient approach. It is especially 
impractical in simulations of stellar collapse, where it is most strategic to restrict the computational domain to just 
outside the stellar surface. CCE is a different approach which is specifically tailored to study radiation at I + . 

In problems with isolated sources, the radiation zone can be compactified inside a finite grid boundary with the 
metric rescaled by 1/r 2 as an implementation of Penrose's [26| conformal boundary at I + . Because I + is a null 
hypersurface, no extraneous outgoing radiation condition or other artificial boundary condition is required. In CCE, 
Cauchy data on the extraction worldtube provides the inner boundary data for a characteristic evolution extending to 
a compactified T + , where the waveform is defined unambiguously by geometric methods. This eliminates waveform 
error due to asymptotic approximations and gauge ambiguities introduced by the choice of extraction worldtube. 
In addition, the extraction worldtube can be placed in the near zone surrounding the sources in order to enhance 
computational efficiency. See [27] for a review. 

High accuracy waveforms from a binary inspiral and merger are important for the design of the detection templates 
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that are critical for the success of gravitational wave astronomy. This has stimulated efforts to increase the accuracy of 
characteristic evolution for use in CCE. Another global approach applicable to isolated systems is to base the Cauchy 
problem itself on the analogue of the hyperboloidal Cauchy hypersurfaces in Minkowski space, which asymptote to 
I + . This approach, first extensively developed by Friedrich 28], is potentially the basis for a very attractive numerical 
approach to simulate gravitational wave production. For reviews of progress on the numerical implementation see 129- 
l31j ] . In spite of the attractiveness of the hyperboloidal approach and its recent success with model problems [32h35| , 
considerable work remains to make it applicable to systems of astrophysical interest. 

The Cauchy evolution codes have incorporated increasingly sophisticated numerical techniques, such as mesh re- 
finement, multi-domain decomposition, pseudo-spectral collocation and high order (in some cases eighth order) finite 
difference approximations. Work has begun to incorporate such techniques in characteristic codes [36|. However, such 
high accuracy methods cannot by themselves cure some of the major complications and sources of error arising in 
CCE. The timelike extraction worldtube T at the inner boundary of CCE is constructed from a coordinate sphere 
x 2 + y 2 + z 2 = 1Z 2 , 1Z = const, cut out from the Cartesian Cauchy grid. However the radial grid points of the 
Bondi-Sachs system are based upon an areal coordinate r, with the angular grid lying on the spheres r = R = const. 
As a result, the extraction worldtube T does not lie on the gridpoints of the Bondi-Sachs system (except for special 
cases such as spherical symmetry) . This necessitates the introduction of an auxiliary characteristic coordinate system 
in the neighborhood of T in which the radial coordinate is replaced by an affine paramter A along the outgoing null 
rays. By taking advantage of the affine freedom, T can then be parametrized by A = 0. The Cauchy data is first 
transformed into the affine characteristic system and expanded about A = to a sufficient power of A to determine 
data for the inner r-grid points of the Bondi-Sachs system in the neighborhood of T . This is a complicated proce- 
dure which introduces interpolation error and has even led to inconsistent inner boundary conditions in the initial 
implementation of CCE (see 10] for a discussion). 

In view of this, the question naturally arises: Why not use the affine-null system in the first place for the char- 
acteristic evolution algorithm and grid? The history behind this choice goes far back. It has to do with the simple 
hierarchical structure that the Einstein equations take in the Bondi-Sachs system, but which is seemingly broken in 
the affine system. We explain this in SecHU In Sec lIIIl we discuss the details of the Einstein equations in the affine-null 
system and show how the hierarchical structure can be restored. This provides the basis for a new worldtube-nullcone 
evolution algorithm. 



II. NULL SPHERICAL COORDINATE SYSTEMS 



The coordinates of both the Bondi-Sachs system and null-affine system are based upon a family of outgoing null 
hypersurfaces emanating from the spherical cross-sections of a timelike worldtube T , where the null coordinate u 
labels these hypersurfaces and the angular coordinates x A (A = 2, 3) label the spherical set of null geodesic rays. In 
the Bondi-Sachs system, the surface area coordinate r labels the points along the outgoing null rays. In the resulting 
x a = {u,r,x A ) coordinates, the metric takes the Bondi-Sachs form (TBI, [ill] 

ds 2 = - (e 213 ^ - r 2 h AB U A U B ^j du 2 - 2e 2f3 dudr - 2r 2 h AB U B dudx A + r 2 h AB dx A dx B , (2.1) 

where det(/i J 4s) = det(q , J 4s) = l{x c ), with qab{x c ) some standard choice of unit round-sphere metric. The fields j3, 
U A , V and Iiab are functions of (u, r,x A ). Here Hab is the metric of the topological 2-spheres (u = const, r = const) 
after conformal rescaling by 1/r 2 to surface area 47r. Its inverse is defined by h AC hcB = S B . 

The affine-null system is similarly based upon the outgoing null hypersurfaces u = const emanating from T with 
coordinates x A labeling the null rays but now an affine parameter A is used to coordinatize points along the rays. The 
affine freedom 

A -> A(u, x A )X + B(u, x A ) 

is used to prescribe the normalization (V a u)V a A = —1 and set A = on T. In the resulting x a = (u, X,x A ) 
coordinates, the metric takes the form 

ds 2 = -(V - g AB W A W B )du 2 - 2dud\ - 2g AB W A dudx B + g AB dx A dx B (2.2) 

In addition, we again set gAB — r 2 hAB, where det Hab = det qAB with qAB(x c ) a unit round-sphere metric. However, 
r is now a metric function of (u, A, x A ) along with W A , V and Hab- 

At the metric level, such affine-null coordinates were introduced by Sachs [371 ] in formulating a double-null initial 
value problem. They are also the natural coordinates adopted in the Newman-Penrose [2fj| formulation of the Einstein 
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equations in terms of a null tetrad and the associated Weyl tensor components. The affine coordinate A only becomes 
singular at caustics whereas the areal coordinate r also becomes singular at points where the expansion of the null 
rays vanish. In particular, this occurs at the points on a stationary event horizon. As a result, codes based upon 
an areal coordinate have poor accuracy in tracking the late time tail preceding black hole formation, cf. [38] for a 
discussion in the context of black hole perturbation theory. The affine null metric and and the Bondi-Sachs metric 
are related by the transformation X(u, r, x A ) determined by 

d r X{u,r,x A ) = e 20 . (2.3) 

However, the simplicity of this transformation is misleading because the surfaces r — const which detemine the 
partial derivatives Da and d u in the Bondi-Sachs system differ from the A = const surfaces that determine the partial 
derivatives 8a and 8 U in the affinc-null system. It is important to keep this distinction in mind in the following 
comparison of the corresponding evolution systems. 

The role of the different components of the Einstein equations in formulating a characteristic initial value prob- 
lem can be best described in terms of an orthonormal null tetrad (L a , N a , M a , M a ), corresponding to the metric 
decomposition 

gab = -L (a N b) + M (a M b) , N a L a = -2, M a M a = 2. (2.4) 

We choose L a = —u a to be the future pointing normal to the null hypersurfaces (so that L a is tangent to the outgoing 
rays) and choose M a to be a complex spatial vector tangent to the null hypersurfaces. This uniquely determines N a . 
Then the vacuum Einstein equations G a b — decompose into the main equations 

L b G ab = (2.5) 
M a M b G ab = (2.6) 



and the supplementary equations 



N b R a b = 0, G a b — Rab — ^9abR- (2.7) 



It is a consequence of the Bianchi identities that if the main equations are satisfied then N R a b satisfies a first 
order ordinary differential equation along the null rays. As a result, if the main equations are satsified and the 
supplementary equations N b R a b = are satisfied on the worduibe T then they will be satisfied everywhere. This 
result was first demonstrated for the Bondi-Sachs system in [TEI, [T^j but it also holds for the affine-null system. 
See [39| for a recent discussion of the supplementary equations as a system of worldtube conservation laws that 
impose symmetric hyperbolic constraints on the worldtube data. In CCE, the worldtube data is supplied by solutions 
of the Einstein equations determined by the Cauchy evolution and it is assumed this data is consistent with the 
supplementary equations. Thus we concentrate here on the main equations. 

First consider the Bondi-Sachs equations. In that case, following the formalism developed in [40(, the main equations 
(|2.5p take the schematic form of hypersurface equations 

/8,r = Mp[hcD] (2-8) 
(^c~ 2f3 h AB U B r )^ = Nu[hcD,P] (2.9) 
V r = Nv[hcD,(3,U c ], (2.10) 



where a "comma" denotes partial derivatives, e.g. /3 r = d r /3, and the main equations (|2.6[) take the form of evolution 
equations 

M A M B (rh AB>u ) >r = N h [h CD , (3, U c , V}. (2.11) 

Here the TV-terms on the right hand sides of (12.81) - (|2.1ip can be calculated from the values of their arguments on 
a given u — const null hypersurface. Moreover, each TV-term only depends upon previous members in the sequential 
order [hcD, ft, U c , V}. Because of this hierarchical structure of the system, given Hab on an initial null hypersurface 
u = 0, the main equations can be integrated radially in sequential order to determine the initial values of /3, U A , V 
and liAB.u at u = in terms of their integration constants on T, i.e. 

(3\ T , U A \ r , U A \ r , V\ T , h AB , u \ T . (2.12) 

In addition, the location of the worldtube, specified by r\f = R(u,x A ), is another essential part of the data. After 
determining Hab,u at u = 0, the hypersurface data Hab can be advanced to u = Au by a finite difference procedure. 
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Given the worldtube data (|2.12[) . this procedure can be iterated to form a worldtube-nullcone evolution algorithm. 
This evolution algorithm is extremely simple and economical compared to Cauchy evolution algorithms. It is the 
algorithm underlying the PITT null code. 

Now consider the afhne-null system, for which the main equations take the schematic form, 



r -1 r,AA - Hr[hcD] (2.13) 

{r 4 h AB W B x ),x = n w [h CD ,r] (2.14) 

2{r 2 ). u - V(r 2 ), A ) = Uv[h C D,r,W c ) (2.15) 

M A M B (rh AB ), uX = H h [h CD ,r,W c ,V], (2.16) 



where the 'H-terms on the right hand sides of (|2.13p - (|2.16[) can again be calculated from the values of their arguments 
on a given u — const null hypersurface. 

As in the Bondi-Sachs case, the "W-terms depends upon the metric functions in sequential order, in this case in the 
order [hcD> r, W , V]. However, the hierarchical structure of the radial integration scheme is broken by the appearance 
of the term (r 2 ) jM term on the left hand side of (|2.150 . Thus (|2.15l) is not a pure hypersurface equation and the radial 
integration scheme does not produce an evolution algorithm in the same was as for the Bondi-Sachs system. This was 
this reason that the afhne-null formulation was not chosen in building the PITT null code. 

However, by reformulating the hypersurface equation (|2.15[) by the introduction of an auxiliary variable, the pure 
hypersurface form of the radial integration scheme can be restored. This new formulation is described in Sec. Mil after 
presenting the details of the main equations of the null-affine system. 

III. THE NULL AFFINE EVOLUTION SYSTEM 

In presenting the details of the Einstein equations for the afhne-null system, we begin with some useful formula for 
describing the metric (|2.2p and its associated connection and curvature. We then proceed to describe the construction 
a numerical evolution algorithm. 

A. Calulation of the Einstein tensor 

The contravariant components of the metric (12. 2[) are given by 

g uX = -l, g uu =g uA = Q, g XA = -W A , g XX =V, g AB =r- 2 h AB . (3.1) 

It is convenient to introduce a dyad vector m A to represent the 2-metric by 

h AB = m (A rh B) , h AB = m ( - A rh B} , m A = h AB m B h AB m A m B = 2. (3.2) 

For that purpose, we chose the vector M a forming the null tetrad (|2.4p to lie tangent to the surfaces (u = const, r = 
const), so that it has components M a = (0,0, M A ). We then set M A — r~ 1 m A . Here we raise and lower indices 
of 2-dimensional vector and tensor fields on the sphere with h AB and h AB , e.g. W A = h AB W B . We recall that 
det h AB = detq AB , where q AB (x ) is some standard choice of unit round-sphere metric. The determinant condition 
implies 

h AB h AB . u = m A rh B h AB , u = 0, h AB h AB . x = m A fh B h AB . x = 0. (3.3) 
We fix the spin rotation freedom in the dyad m A — > e"p( u <*> x )m A by requiring 

fh A m A \ = (3.4) 

and 

m A m A , u \x =0 = 0. (3.5) 
Note that the determinant condition (|3.3[) also implies 

m m At x + m A fh A ^\ — 0. (3-6) 
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The spin rotation freedom then reduces to a phase factor <p(x A ), which is determined by the choice of conventions at 
(u = 0, A = 0). Given these conventions, h A B and m A are in one-to-one correspondence. 
As an example, for stereographic coordinates x A = (77, p) on the unit round-sphere with metric 



q AB dx A dx B = ^q(dr} 2 + dp 2 ), y/q = — - T — (3.7) 

1 + if + P 

the rescaled metric on the general curved topological sphere can be represented as 

, _ r- ( e 27 cosh 2a sinh2a \ , , 

h *B-Vq{ sinh2a e -27 COsh2a J' (3 - 8) 

iAB 1 /" e -27 cosh 2a — sinh2a \ , . 

~ ^/f ^ -sinh2a e 27 cosh 2a J - ^ j 

where 7 and a represent the two degrees of freedom. A specific choice of polarization dyad associated with this 
representation is 

m A = q 1 ^ ^e 7 (cosh a + i sinha), ie _7 (cosha — i sinha)^ , (3.10) 

m A = q^ 1 / 4 I e~ 7 (cosha — i sinha), ie 7 (cosha + i sinha) I . (3-11) 



The components of the Einstein tensor can be calculated in terms of the metric functions h A B, t, W a , V from the 
components of the curvature tensor 

R a bcd = d c T a bd — ddT a bc + r a ce r e M — r a ci e r e f, c , R a b — R C a cb, (3.12) 

where 

n c = \g ad {d b9cd + d c9bd - d d9bc ) (3.13) 

are the Christoffel symbols. The components of the Christoffel symbols in terms of the metric functions are given in 
Appendix [A] We denote 2-dimensional covariant derivatives of tensor fields on the sphere with respect to Hab by a 
"colon", e.g. 

W A -b = d B W A + ^T A C W C (3.14) 

where 

(h) ric = \h AD {d B h CD + d c h BD - d D h BC ) (3.15) 

is the Christoffel symbol associated with h A B- 

In terms of these conventions and notation, the main affine-null equations (|2.13l) - (12.16[) have the specific form 

= R{u b = ^ - \h c x D h C D,x (3.16) 
r 4 ' 

= R\u h = -^-{^h AB W B x ), x + (^), A - ( r -^)h BC h AC ,x - \h BC h AC ,X:B (3.17) 
2r £ ' r r 2 

= h AB R AB = {2) n + (V 2 X„ - V(r 2 ), A ) + (±(r*W A )^ - (lnr 2 y A :A - r \ AB W A W B (3.18) 
= m A m B R AB = m A m B {r(rh AB ), u x - ^(r 2 Vh AB ,x)^x + rr :C W c h AB ,x + (r 2 W A:B ),x 

2 2 4 

+ ^h AB ,\m c m D (W C:D - W D:C ) + r —W c h AB ,x:C ~ ^h AC h BD W c x W^}. (3.19) 
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B. Restoration of the hypersurface equation hierarchy 

The strategy now is to use the auxiliary variable 

y = V-^ (3.20) 

to eliminate the explicit appearance of thr r u derivative in (13.25)) and reexpress it as a hypersurface equation for y. 
In that process, substitution of y for V in the evolution equation (13. 19)) leads to the intermediate expression 

m A m B {r(rhAB),u\ - \{r 2 VhAB,x) .a} = rm A m E ''{rh,AB,u - — ^/MlB,a},A ~ m A m B \{r 2 yhAB,x),x (3-21) 
2 r, a 2 

or, using (|3.4[) . 

1 r 1 

m A m B {r(rhAB),u\ ~ 7;(r 2 VhAB,x).x} = r{m A m B r(h A B,u — hAB,x)}.x-m A m B -{r 2 yh A B,x),x- (3.22) 

2 r,A 2 

When (13.221) is inserted back into (|3.19p it gives a hypersurface equation for a combination of the time derivatives 
m A m B Hab,u and r u . In order to complete an evolution system, we need an additional radial equation for r „. This 
is obtained from the w-derivative of the Raychauduri equation (|3.16[) . which determines the rate of change of the 
expansion of the outgoing rays. In order to formulate a hierarchical radial integration scheme, we introduce the 
auxiliary variables 



and 



Then the w-derivative of ()3.16p gives 



P = r, u . (3.23) 

kAB = h A B,u- (3.24) 

P,xx = ^h% D h CD ,x + ^h B x c k BC ,x, (3.25) 

where the determinant condition implies that the undfferentiated kAB terms vanish, i.e. 

h BD h CB h ED ^k BC = 0. (3.26) 

In assembling the foregoing results into a radial integration hierarchy, simplifications result from using the spin- 
weighted scalars 

a = ^m A m B h A B,\ = ^m A m A ,x , (3.27) 

K = — m A m B kAB, kAB — nrhAifiB + srajniB (3.28) 



and 



Then (pT2"2"j) becomes 



J = 4(r a k - per) = r,\m A m B h A B,u - r iU m A m B h A B.x- (3.29) 



m A m B {r{rhAB),uX-Ur 2 Vh AB ,x),x} = r(— ) - m A m B hr 2 yh A B,x),x (3.30) 

and p. 251) becomes 

P,XX = -paa - raK,x - raR t x 

= -^(raa),x ~ 2raa( — )\ - ^{-),x - ^(-),a (3.31) 
r,x r,x 4 r,\ 4 r,\ 

where, from fXT^) and (J3H7J), 

raa = -r.AA- (3.32) 
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C. The evolution algorithm 

By assembling the results in Sec. IIIIB1 the main equations f!3. 16[) - (|3.19l) . along with p.31[) . now take the desired 
form 

r _1 r,AA = -oo = H r [hcD\ (3.33) 

(r 4 h AB W B x ). x = 2r 2 ( r -^). A -2rr, B h BC h AC , x -r 2 h BC h Ac .X:B 
r 

= H w [h CD ,r] (3.34) 
3VXa) = m n+(±(r*W A ),>?) -(inr 2 )^- r ^h AB Wf x W B x 

= fly|7»ciJ,r,W c ] (3.35) 

— ) - m A m s { T i(r 2 ^ AB , A ), A -r :C iy c /i AB , A -i(r 2 ^ :B ), A 
rW > 2r r 



^h AB ^rh c m D (W c - D - IU D:C ) - ^ c ^s,A:C + \h Ac h BD W c x W D x } 
Hj[h C D,r,W c ,y] (3.36) 
/M ' (a{—) >x + a{—) >x \ =H p [h CD ,r,J], (3.37) 



r ,A/ , AA 4?\ A V r A ' r x 

where the H -terms can be calculated from the values of their arguments on a u = const null hypersurface. Given 
h A B {u, A, x ), the system (|3 .33[) - (|3 .37[) forms a hierarchy of radial hypersurface equations which can be integrated 
to determine the remaining variables in the sequential order [r, W A , y, J, p]. 

This gives rise to the following initial-boundary value problem. Specify the initial hypersurface data 

[h AB ,r], u = 0, A>0 (3.38) 

and the initial boundary data 

[r,x,W A ,W£,y,J t p,p,x], u = 0, A = 0, (3.39) 

subject to the constraint (|3.33l) . which constitutes an ordinary differential radial equation. On the boundary, specify 

[W A ,W^,y,J,p,p, X ], u>0 : A = 0. (3.40) 

Using the initial data, integrate (|3.33[) - (|3.37[) to determine the initial values of [r, W A ,y, J, p\. 

Given this initialization and the boundary data, the evolution system can be integrated by a finite difference 
approximation. The initial values p(0,X,x A ) and J(0, X 7 x A ) determine the values of r and h AB at u = Au through 
p.23|) and (|3. 291) . Using the boundary data, (|3. 341) - p. 371) can then be integrated in sequential order to determine 
[W A ,y,J,p] at«= Au. Now [h AB ,r,W A ,y,J,p] are known at u = Au and this process can be repeated to provide 
a finite difference evolution algorithm. If the algorithm converges as Am — > then it produces a solution to the 
affine-null initial-boundary value problem for Einstein's equations. 



IV. DISCUSSION 



We have constructed an evolution algorithm based upon the affine-null system which, like the Bondi-Sachs system, 
is based upon a hierarchy of radial equations along the outgoing characteristics. It has the additional advantage of the 
flexibility in describing an arbitrary inner worldtube boundary as a coordinate surface. This is especially important for 
application to CCE, where the inner boundary, which is constructed in terms of the Cauchy coordinates, generically 
differs from the r = const Bondi-Sachs worldtubes. As a result, the affine-null algorithm offers the possibility of 
increased economy and accuracy. 

The formal solution of the null-affine problem constructed in Sec. IIII CI yields an exact solution provided the finite 
difference approximation converges as Au — > 0. A necessary condition for this is the well-posedness of the underlying 
analytic initial-boundary value problem. Well-posedness, i.e. the existence of a unique solution which depends 
continuously on the data, is a necessary condition for a successful numerical treatment. Although characteristic 
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evolution codes based upon the Bondi-Sachs formalism have been demonstrated to be stable in a large number of 
test cases [H, Hl| , there remains some lingering doubt because well-posedness of the analytic problem has not yet 
been established. Rendall [42[ has shown that the affine-null problem is well-posed in the double null case where 
the inner boundary is also a null hypersurface. However, Rendall's approach cannot be applied to the corresponding 
problem where the inner boundary is a timelike worldtube. The well-posedness of the worldtube-nullcone characteristic 
initial-boundary value problem for Einstein's equations remains an outstanding issue. 

The only source of error in CCE which does not decrease with numerical resolution arises from the mismatch 
between the initial Cauchy data and initial characteristic data. This results because the radius 72.2 of the outer Cauchy 
boundary is larger than the radius IZi of the inner worldtube boundary of the charcteristic evolution. Whereas the 
Cauchy data in the region 7?.! < 1Z < H2 is chosen, say, by some constraint solver for binary black hole initial data, 
the characteristic initial data for 1Z > IZi is chosen to suppress the initial radiation content by requiring that the 
Newman-Penrose Weyl component ^0 = 0. This mismatch between the initial Cauchy and characteristic data leads 
to an extraneous error in the extracted waveform which is related to the spurious radiation content in the Cauchy 
data. This error decreases as IZ2 — > Hi but present day Cauchy codes require that the outer boundary be in the far 
zone of a binary black hole to avoid incoming radiation generated by the outer boundary condition. As a result, the 
full potential of CCE is not realized. 

This mismatch can be eliminated by Cauchy-characteristic matching (CCM) (43[. CCE is one of the pieces of CCM 
in which the characteristic worldtube data is extracted from the Cauchy evolution. In CCM, data on the outer Cauchy 
boundary is in turn obtained from the characteristic evolution. In doing so, it is possible to place the radius IZ2 of the 
Cauchy boundary just outside the radius 1Z\ of the characteristic extraction worldtube. In fact, in a finite difference 



implementation of CCM for a model scalar wave problem [44|, it has been possible to arrange that IZ2 — > Hi in 
the continuum limit. This resulted in a seamless interface between the Cauchy and characteristic evolutions with no 
mismatch in the initial data. 

The success of CCM depends upon the proper mathematical and computational treatment of the initial-boundary 
value problem (IBVP) for the Cauchy evolution. At present, the only successful 3D application of CCM in general 
relativity has been to the linearized matching problem between a characteristic code and a Cauchy code based upon 
harmonic coordinates [45j . Considerable work remains to apply it to astrophysical systems. The linearized harmonic 
code satisfied a well-posed initial-boundary value problem, which seems to be a critical missing ingredient in earlier 
attempts at CCM in general relativity. More recently, a well-posed initial-boundary value problem has been established 
for fully nonlinear harmonic evolution [46l. |47] . which should facilitate the extension of CCM to the nonlinear case. 
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Appendix A: Christoffel symbols 

The calculation of the Ricci tensor 

Ra/3 — ^a/3,u — ^va,P + ^afi^'pv ~~ ^uoFp0 (Al) 

which enter the main equations can be carried out explicitly in terms of the Christoffel symbols for the affine-null 
metric (|2.2p . The components are listed according to the notation (x°,x 1 ,a; j4 ) = (u, X,x A ). 

T a al/ = d u \n{r 2 ^q) (A2) 
r? a = (A3) 
r?i - (A4) 



T AB = rr ,\ h AB + -^h AB .x (A5) 
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^ab = -jihAcWg + h BC W c A + h AB . u - Vh AB , x ) + rh AB (r tC W c + r „ - Vr A ) (A6) 

rSs = rr, xW c h AB + \r 2 W c h AB + + -(r, B 8 c A + r. A 5 c B - r, D h CD h AB ) + ^T C AB (A7) 

Z V 

K A = -^{h AB . x W B + h AB W B x ) - rr. x h AB W B (A8) 

rio = \v,x - ^h AB W A W B x (A9) 

= ^h AB W B x (A10) 

^ = -\v,x + \(r 2 h AB W A W B ), x (All) 

rSi = ~r,x5 c A + V'Wa (A12) 
r z 

rj = i^, + \{r 2 h AB ) fi W A W B + iv(V - r 2 h AB W A W B ), x - \w A (V - r 2 h BC W B W c ). A (A13) 
Tqa = + \v(r 2 h AB W B ), x - rr. B W B h AC W c - y W^AbcWS + h AC Wg) - ^W B (r 2 h AB ) fi (A14) 

rg* - - V B (r 2 ^c^ C ),A + ^ BC (r 2 /uc),o - — W B - + I^/^ AC T4^ + (A15) 

Z Z T Z T Z 

T a = _ r _A W A _ l h AB (h BC W c ). x (A16) 

V z 

T A = - l -W A V++ l -W A {r 2 h BC W B W c ), x -2W Ar -f - h AB (h BC W c ), u + ^h AB V B 

r -^h AB h CD W c W D -\h AB (h CD W c W D ), B . (A17) 



[1] "Cauchy-characteristic extraction in numerical relativity", N.T. Bishop, R. Gomez, L. Lehner, and J. Winicour, Phys. Rev. 
D 54 6153 (1996). 

[2] "The NINJA-2 catalog of hybrid post-onian/numerical-relativity waveforms for non-precessing black-hole binaries", P. 

Ajith et al: NINJA project, Class. Quantum Grav. 29 124001 (2012). 
[3] "Unambiguous determination of gravitational waveforms from binary black hole mergers", C. Reisswig, N. T. Bishop, D. 

Pollney and B. Szilagyi, Phys. Rev. Lett. 103, 221101 (2009). 



10 



[4] "Characteristic extraction in numerical relativity: binary black hole merger waveforms at null infinity", C. Reisswig, N. T. 

Bishop, D. Pollney and B. Szilagyi, Class. Quant. Grav. 27, 075014 (2010). 
[5] "Binary black hole waveform extraction at null infinity", M. C. Babiuc, J. Winicour and Y. Zlochower, Class. Quant. Grav. 

28, 134006 (2011). 

[6] "Gravitational wave extraction in simulations of rotating stellar core collapse", C. Reisswig, C. D. Ott, U. Sperhake and 

E. Schnetter, Phys. Rev. D 83, 064008 (2011). 
[7] "Dynamics and gravitational wave signature of collapsar formation", C. D. Ott, C. Reisswig, E. Schnetter, E. O'Connor, 

U. Sperhake. F. LofHer, P. Diener, E. Abdikamalov, I. Hawke and A. Burrows, Phys. Rev. Lett. 106, 161103 (2011). 
[8] "Gravitational memory in binary black hole mergers", C. Reisswig and D. Pollney, Astrophys. J. Lett. 732, L13 (2011). 
[9] "Gravitational- wave detectability of equal-mass black-hole binaries with aligned spins", C. Reisswig, S. Husa, L. Rezzolla, 

E. N. Dorband, D. Pollney and J. Seiler, Phys. Rev. D 80, 124026 (2009). 
[10] "A characteristic extraction tool for gravitational waveforms", M. C. Babiuc, B. Szilagyi, J. Winicour and Y. Zlochower., 

Phys. Rev. D 84 044057 (2011). 
[11] http:/ /www. einsteintoolkit.org 

[12] "Null cone computation of gravitational radiation", R.A. Isaacson, J.S. Welling and J. Winicour, J. Math. Phys. 24 1824 
(1983). 

[13] "High-powered gravitational news", N. T. Bishop, R. Gomez, L. Lehner, M. Maharaj and J. Winicour, Phys. Rev. D 56 
6298 (1997). 

[14] "Gravitational fields in finite and conformal Bondi frames", L. A. Tamburino and J. Winicour, Phys. Rev. 150 1039, 1966. 
[15] "Gravitational waves in general relativity VII. Waves from axi-symmetric isolated systems", H. Bondi, M. van der Burg 

and A. W. K. Metzner, Proc. R. Soc. London A 269, 21 (1962). 
[16] "Gravitational waves in general relativity VIII. Waves in asymptotically flat space-time"', R. K. Sachs, Proc. R. Soc. 

London A 270, 103 (1962). 
[17] "Stability of a Schwarzschild Singularity", T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 
[18] "Gravitational field of a particle falling in a Schwarzschild geometry analyzed in tensor harmonics", F. Zerilli, Phys. Rev. 

D. 2, 2141 (1970). 

[19] "Gauge-invariant non-spherical metric perturbations of Schwarzschild black-hole spacetimes", A. Nagar and L. Rezzolla, 

Class. Quantum Grav. 22, R167 (2005). 
[20] "An approach to gravitational radiation by a method of spin coefficients", E. T. Newman and R. Penrose, J Math. Phys. 

3, 566 (1992). 

[21] "Modeling gravitational radiation from coalescing binary black holes", J. Baker, M. Campanelli, CO. Lousto and R. 

Takahashi, Phys. Rev. D, 65, 124012-124034 (2002). 
[22] "Evolution of binary black-hole spacetimes", F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005). 

[23] "Accurate evolutions of orbiting black-hole binaries without excision", M. Campanelli, C. O. Lousto, P. Marronetti and 

Y. Zlochower, Phys. Rev. Lett. 96 111101 (2006). 
[24] "Binary black hole merger dynamics and waveforms", J. G. Baker, J. Centrella, D-I. Choi, M. Koppitz and J. van Meter, 

Phys. Rev. D 73, 104002 (2006). 
[25] "Dealing with delicate issues in waveform calculations", L. Lehner and O. M. Moreschi, Phys. Rev. D 12, 124040 (2007). 
[26] "Asymptotic properties of fields and space-times", R. Penrose, Phys. Rev. Lett. 10, 66 (1963). 
[27] "Characteristic evolution and matching", J. Winicour, Living Rev. Rel. 15, 2 (2012). 

[28] "Cauchy problems for the conformal vacuum field equations in general relativity", H. Friedrich, Commun. Math. Phys. 91, 
445 (1983). 

[29] "Conformal infinity", J. Frauendiener, Living Rev. Rel. 7, 1 (2004). 

[30] "Problems and successes in the numerical approach to the conformal field equations", S. Husa, Lecture Notes Phys. 604, 
239 (2002). 

[31] "Numerical relativity with the conformal field equations", S. Husa, Lecture Notes Phys. 617, 159 (2003). 
[32] "Hyperboloidal evolution with the Einstein equations", A. Zenginoglu, Class. Quantum Grav. 25, 195025 (2008). 
[33] "Regularity of the Einstein equations at future null infinity", V. Moncrief and O. Rinne, Class. Quantum Grav. 26, 125010 
(2009). 

[34] "An axisymmetric evolution code for the Einstein equations on hyperboloidal slices", O. Rinne, Class. Quantum Grav. 27, 
035014 (2010). 

[35] "Tetrad formalism for numerical relativity on conformally compactified constant mean curvature hypersurfaces" , J. M. 

Bardeen, O. Sarbach and L. T. Buchman, Phys. Rev. D 83, 104045 (2011). 
[36] "General relativistic null-cone evolutions with a high-order scheme" C. Reisswig, N. T. Bishop, and D, Pollney, 

larXiv:1208.389TI (2012). 

[37] "On the Characteristic Initial Value Problem in Gravitational Theory", R. K. Sachs, J. Math. Phys. 3, 908 (1962). 

[38] "Retarded radiation from colliding black holes in the close limit" S. Husa, Y. Zlochower, R. Gomez and J. Winicour, 

Phys. Rev. D 65 084034 (2002). 
[39] "Worldtube conservation laws for the null-timelike evolution problem", J. Winicour, Gen. Rel. Grav. 43, 3269 (2011). 
[40] "Newtonian gravity on the null cone" J. Winicour. J. Math. Phys. 24, 1193 (1983). 

[41] "Moving black holes in 3D", R. Gomez, L. Lehner, R. Marsa and J. Winicour, Phys. Rev. D 57, 4778 (1998). 

[42] "Reduction of the characteristic initial value problem to the Cauchy problem and Its applications to the Einstein equations" , 

A. D. Rendall, Proc. R. Soc. London A, 427, 221 (1990). 
[43] "Cauchy-Characteristic Matching", N.T. Bishop, R. Gomez, L. Lehner, B. Szilagyi, J. Winicour and R. A. Isaacson, in B. 



11 



Iyer and B. Bhawal (Eds.), Black Holes, Gravitational Radiation and the Universe, Kluwer Academic Publishers, Dordrecht 
(1998). 

[44] "Cauchy-characteristic evolution and waveforms", N. T. Bishop, R. Gomez, P. R. Holvorcem, R. A. Matzner, P. Pa- 

padopoulos and J. Winicour, J. Comput. Phys. 136, 140 (1997). 
[45] "Well-Posed Initial-Boundary Evolution in General Relativity" B. Szilagyi and J. Winicour Phys. Rev. D 68, 041501 (2003). 
[46] "Problems which are well-posed in a generalized sense with applications to the Einstein equations", H.O. Kreiss and J. 

Winicour, Class. Quantum Grav. 23, S405-S420 (2006). 
[47] "Boundary conditions for coupled quasilinear wave equations with application to isolated systems", H-O. Kreiss, O. Reula, 

O. Sarbach and J. Winicour, Commun. Math. Phys. 289, 1099 (2009). 



